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A general numerical solution of the dispersion integral relation between the real and the imag- 
inary parts of the nuclear optical potential is presented. Fast convergence is achieved by means 
of the Gauss-Legendre integration method, which offers accuracy, easiness of implementation and 
generality for dispersive optical model calculations. The use of this numerical integration method 
in the optical-model parameter search codes allows for a fast and accurate dispersive analysis. 
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For many years the evaluation of reaction cross section and elastic scattering data has relied on the use of the 
optical model. A significant contribution during the last two decades can be considered the work of Mahaux and 
co-workers on dispersive optical-model analysis |I] 0]. The unified description of nuclear mean field in dispersive 
optical- model is accomplished by using a dispersion relation (DR) , which links the real and absorptive terms of the 
optical model potential (OMP). The dispersive optical model provides a natural extension of the optical model derived 
data into the bound state region. In this way a physically self-consistent description of the energy dependence of the 
OMP is obtained and the prediction of single-particle, bound state quantities using the same potential at negative 
energies became possible. Moreover additional constraint imposed by DR helps to reduce the ambiguities in deriving 
phenomenological OMP parameters from the experimental data. The dispersive term of the potential AV(E) can be 
written in a form which is stable under numerical treatment [Q , namely: 



AV(F) - 2 (F F ) f°° W(E')-W(E) 
AV(E) ~ * iE ~ Ef) J Ef (E' — Ep) 2 ~ (E — Ep) 



(1) 



where E is the incident projectile energy and Ep is the Fermi energy for the target system. 

The imaginary OMP can be written as the sum of volume and surface contributions: W(E) — Wy(E) + Ws(E). 
It is useful to represent the variation of surface Ws(E) and volume absorption potential Wv{E) depth with energy 
in functional forms suitable for the DR optical model analysis. An energy dependence for the imaginary volume term 
has been suggested in studies of nuclear matter theory: 

where Ay and By are undetermined constants. Mahaux and Sartor 0] have suggested n = 4, while Brown and Rho 
[?], n = 2. An energy dependence for the imaginary-surface term has been suggested by Delaroche et al 0] to be: 

W S (E) = A s (E -E F ) m ,_ c > E _ S h (3) 

*v ) (E - E F ) m + {B s ) m K 1 u w 

where m — 4 and As, Bs and Cs are undetermined constants. Other energy functionals are suggested in the literature 
According to equations @ and (||) the imaginary part of the OMP is assumed to be zero at E = Ep and nonzero 
everywhere else. A more realistic parametrization of Wy(E) and Ws(E) forces these terms to be zero in some region 
around the Fermi energy. A physically reasonable energy for defining such a region is the average energy of the 
single-particle states Ep We can define the offset-energy E a f for both the particle and the hole region by the 
relation: 

E of = E P -E F (4) 
Therefore a new definition for imaginary part of the OMP can be written as: 

Wy(E) = 0, for E F <E<E p (5) 

, . Ay(E-E P ) n 

Wv(E) = t= — — r ' , , for Ep < E 6 

v ; (E - E P ) n + (B v ) n y ' 



1 



and likewise for surface absorption. The symmetry condition W(2Ep — E) — W(E) is used to define imaginary part 
of the OMP for energies below the Fermi energy. We assume below that equations (g) and (§) are modified according 
to (|) and (§). 

In a recent work || , analytical solutions of the dispersion integral relation between the real and imaginary parts of 
the nuclear optical model were proposed. The main reason for such effort was the slow convergence of the Simpson 
integration method (SIM) for solving the DR integral (0). The authors needed 6000 steps in total (50 keV each) to 
calculate (Q) up to 300 MeV. Beyond 300 MeV, the integral was approximated by assuming that the numerator of 
its integrand become constant- valued ||. Obviously this lengthen very much the OMP calculations, so they derived 
approximate analytical solutions of the integral (Q) for usual forms of the imaginary OMP, including those given by (^) 
and (||). However the approximate analytical solution recently derived j(| depends, of course, on assumed functional 
form of the imaginary part of the OMP. 

In this work we propose a fast and general numerical solution for calculation of the DR integral (|l|) , based on Gauss- 
Legendre integration method (GLIM) Q|. Integral in a finite interval (a, b) can be converted by linear transformation 
to the integral in the interval (—1, 1) and then approximated according to the following formula: 



L N 

f{x)dx = Y, WfKX?) (7) 
, * — »&= i 

where WF and Xf are weights and abscissas corresponding to the N point Gauss-Legendre integration method M . 
In this work we were using N=10 for each interval. At the first look of equation ([!]) we could be tempted to use 
Gauss-Laguerre integration method 0. However beyond some energy cut-off E cut , such that E cut >> E, the integral 
can be very well approximated by assuming that the numerator of its integrand becomes constant- valued. In this case 
we obtain the following relation: 

av (e)- 2 (f f)[°° W{E ' ] - W{E) jf> W{E) lor f Ecut - {E - Ef) \ m 

AV res (E) - ~(E - E F )J E ^ {E , Ep)2 _ {e _ ^dE - — log - {E Ep) ) (8) 

Therefore we can apply GLIM from Ep up to E cut for a finite interval. This option was shown to be more accurate 
for a given problem. This method is suitable for practically any type of DR integrand function. Moreover the 
implementation of this method in optical model parameter search code is straightforward. The method is compared 
with SIM and, for the surface dispersive contribution, with the results of analytical integration B. The parameters 
of the dispersive OMP used in our tests correspond to the "partially constrained set" derived in Ref. || to describe 
neutron scattering on 209 Bi up to 80 MeV incident energy (see Table II, ||). For 209 Bi the Fermi energy is equal to 
-5.98 MeV and the offset energy was taken as 3.26 MeV, assuming the average energy of the particle states is the 
same of 208 Pb §. 

In order to attain a good accuracy in the numerical integration (|l|), attention must be paid to the behavior of the 
integrand, 

Iee (E')= W ( E ')- W W (9) 
E ' Ef [ ' (E> - E F Y -{E- E F f 1 ' 

The integrand Ie,e f {E') for dispersive surface correction (Ws(E) given by equation (^) ) is shown in Figure 1. As 
a result of the exponential damping in the surface OMP the integrand Ie,e f {E') decreases sharply above 5-6 MeV. 
This behavior suggests that separation of the integration interval is desirable in order to improve the accuracy of 
the numerical integration method. It should be stressed that exponential damping contained in functional form (|^) 
reflects the physical fact that surface absorption is expected to decrease fast as the incident energy is increased. This 
behavior is observed in almost any phenomenological OMP analysis, so dividing integration interval is not a recipe 
for functional form (||), but should be valid practically for any imaginary surface absorption energy functional. We 
decided to divide the interval in two parts from Ep up to Ei nt = lOMeV and then from Ei nt up to some cut-off 
energy in the integral. The Ei n t position was chosen in order to guarantee monotonous behavior of the integrand 
in the interval from Ei nt up to infinity. Practically it is enough to select Ei nt slightly larger than the energy of the 
maximum value of the integrand. The cut-off energy E cut in the DR integral for dispersive surface correction was 
taken at 1000 MeV. Beyond 1000 MeV we used equation (H) to estimate contribution from the tail of the integral. 
This residual value was added to both numerical methods. The dispersive surface correction AVs(E) calculated by 
GLIM (using 20 integration points) and SIM (100 keV step, more than 10000 integration points) were compared with 
analytical solution from Ref. || . Maximum deviation between any of the three methods was less than 1 KeV in the 
whole energy range from Ep up to 150 MeV energy. 

A similar behavior of the integrand Ie,e f (E') is observed for the case of the dispersive volume correction (Wv(E) 
given by equation (0) ). However the integrand Ie,e f (E') is smoother and has longer tail than was the case for 
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surface correction. The reason is that no exponential damping is introduced here. Similarly to what we did before, 
we divided the interval in two parts from Ep up to Ei nt = 4QMeV and then from Ei nt up to cut-off energy E cut . 
Beyond 1000 MeV we used equation (||) to estimate contribution from the tail of the integral. The integration was 
performed by GLIM and SIM. An excellent agreement is achieved between both numerical methods. 

In conclusion, we have proposed an accurate and general numerical method for calculating dispersive contribution 
to the real part of the optical model potential. The method is very fast and thanks to its computational simplicity 
should be easy to implement in current generation of the optical model parameter search codes. We stress that, unlike 
analytical solution, any energy functional representation of the imaginary part of the OMP can be treated in our 
approach. These properties make the method very useful for introduction of the dispersive optical model relation in 
large scale nuclear data calculations. 
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FIG. 1. Energy dependence of the DR integrand Ie,e f (E') for the dispersive surface contribution of the imaginary OMP 
(short dottted line) E — Ef + 0.25MeV, (dash dotted line) E = 0, (dotted line) E = WMeV 
(dashed line) E = SOMeV, (solid line) E = 50MeV 
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